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We study numerically the influence of disorder and localization effects on the local spectroscopic 
characteristics and infrared optical properties of Gai-^Mn^As. We treat the band structure and 
disorder effects at an equal level by using exact diagonalization supercell simulation method. This 
method accurately describes the low doping limit and gives a clear picture of the transition to 
higher dopings, which captures the localization effects inaccessible to other theoretical methods 
commonly used. Our simulations capture the rich mid-gap localized states observed in scanning 
tunneling microscopy studies and reproduce the observed features of the infrared optical absorption 
experiments. We show clear evidence of a disordered valence band model for metallic samples in 
which (i) there is no impurity band detached from the valence band, (ii) the disorder tends to 
localize and pull states near the top of the valence band into the gap region, and (iii) the Fermi 
energy is located deep in the delocalized region away from the mobility edge. We identify localized 
states deep in the gap region by visualizing the probability distribution of the quasiparticles and 
connecting it to their respective participation ratios. The analysis of the infrared-optical absorption 
data indicates that it does not have a direct relation to the nature of the states at the Fermi energy. 


I. INTRODUCTION 

Gai-^Mn^As is a material system prototype, which 
incorporates carrier-mediated ferromagnetism into 
semiconductorsP® Under equilibrium growth condi¬ 
tions, the solubility of Mn in GaAs is limited to ~ 0.1%. 
Through non-equilibrium low-temperature molecular- 
beam-epitaxy (LT-MBE) technique, high Mn-doped 
samples with more than 1% Mn can be obtained. 6 7 As 
has been shown in a recent systematic study,high 
quality homogenous samples with reproducible charac¬ 
teristics can be prepared by introducing optimization 
of growth and post-growth annealing procedures for 
each doping concentration. Although its potential for 
applications has been curtailed due to growth limitations 
to achieve magnetic transition at room temperature, 
Gai-^Mn^As is a promising material for testing known 
magneto-transport, magneto-optical mechanisms and 
discovering of new magnetic phenomena:^® 

However, even after many recent studies on this mate¬ 
rial, a common debate still exists regarding the electronic 
states near the top of the valence band and at the Fermi 
energy.^ At small Mn doping, x < 0.1%, Gai-^Mn^As 
is insulating and paramagnetic. The bound hole states 
introduced by Mn form an impurity band (IB) detached 
from the top of the GaAs valance band (VB). For higher 
Mn doping, 0.5% < x < 1.5%, the isolated bound hole 
states start overlapping with each other and the IB starts 
mixing with the VB. The material is still insulating but 
ferromagnetism begins to occur. At x ~ 1.5%, the 
abrupt increase of low-temperature conductivity shows 
that the material becomes a degenerate semiconductor. 
In this metallic regime, one proposal states that the IB 
merges with the VB forming a disordered VBP^ An¬ 


other proposal assumes the Fermi level still resides in the 
narrow IB detached from the VBF 7 However, the latter 
proposal, based on early spectroscopic studies on unop¬ 
timized samples,® 7 1 is inconsistent with recent spectro¬ 
scopic studiesPlt has also been shown to be inconsistent 
with microscopic band structure theories? 11 In this paper, 
through a more realistic treatment of disorder effects ex¬ 
trapolated from the small to large doping regimes, we 
show clear evidence of a disordered valence band model 
for metallic samples in which (i) there is no impurity 
band detached from the valence band, (ii) the disorder 
tends to localize and pull states near the top of the va¬ 
lence band into the gap region, and (iii) the Fermi energy 
is located deep in the delocalized region away from the 
mobility edge. 

Some early numerical studies either ignored disor¬ 
der effects or treated them in the framework of mean 
field theory! 1 ®^ In order to treat these effects more 
realistically, we use the exact diagonalization supercell 
method. 18 ^ The advantage of this method is that the 
band structure and disorder effects are captured at an 
equal level which allows the study of disorder and local¬ 
ization effects on the local spectroscopic characteristics 
and infrared optical properties of Gai-^Mn^As. The 
Hamiltonian includes the k • p description of the GaAs 
valance band, the Coulomb interactions at the Hartree 
level, short-range central-cell potential and kinetic ex¬ 
change interaction. By diagonalizing this Hamiltonian 
numerically in the framework of the envelope function 
approximation,^ eigenvalues and eigenfunctions can be 
obtained. Our wave function probability distribution vi¬ 
sualization shows clear hydrogenic-like bound states for 
nearly isolated Mn impurities in low doping samples as 
well as the transition of these bound states to higher dop- 
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ing samples. The (local) density of states (LDoS/DoS) 
calculations capture the rich mid-gap localized states ob¬ 
served in scanning tunneling microscopy studies.^ 2 The 
optical conductivity calculations reproduce the observed 
broad peak at ^ 200 meV of the infrared optical ab¬ 
sorption measurementsAn analysis of the optical ab¬ 
sorption data show no direct relation of this data to the 
nature of the states at the Fermi energy, contradicting 
the implied connection assumed by the IB models. 

The paper is organized as follows. Section II comprises 
a brief introduction of our simulation Hamiltonian. Sec¬ 
tion III provides the simulation results and discussions. 
In this section, we study the bound states properties in 
low Mn-doping limit first, then DoS and LDoS in the 
high concentration regime, followed by the analysis of 
the localized properties of the states, and finally the ac- 
conductivity calculations. Section IV presents the con¬ 
clusions. 


local moment which are aligned in the z direction at zero 
temperature. The phenomenological parameters we use 
are the same as the parameters in Ref. 18] We also treat 
the mutual interaction between holes by finding the self- 
consistent solution of the Hartree potential. We use 6 x 
6x6 nm cube throughout our simulation. 


II. MODEL HAMILTONIAN 


We use a phenomenological model employed by Yang 
et a/. 18 The host semiconductor GaAs valance band is 
described by the six-band k • p Kohn-Luttinger model! 20 
Within the single particle approximation, the hole car¬ 
riers interact with randomly placed Mn local moments 
via the Coulomb and the exchange interaction. A short- 
range central-cell correction is also considered in order 
to capture the difference in the electronegativity of the 
Mn and host Ga atom. 21 In real samples, there are charge 
and magnetic moment compensations due to As anti-sites 
and Mn interstitials. In this paper, we consider only the 
charge compensation introduced by As anti-sites because 
the magnetic moment compensation can be reduced by 
annealing procedures. The total Hamiltonian is given by: 


N M n 


H = H l + Si ■ sJ(r — R,) 


1=1 


N M n 


£ - 

1=1 \ 

N As 

E 


e\r — 
2e 2 


K=1 ^ 1 ^* Ras,K 


- / r o j j 

l (1) 


where H L describes the host valance band, Vo term is the 
central-cell correction, J(f) = J^/[(27rao) 3//2 ]e _r2 / 2a o, / 
is a 6 x 6 unit matrix, and s = (s x ,s y ,s z ) where s X:V:Z 
are 6x6 matrices which describe hole spinsf 22 lRMn,/ and 
Ras,k are the positions of Mn and As respectively. The 
number of holes is given by relation Nh = Nmu ~ 2 Na s • 
Our numerical method diagonalizes the single particle 
Hamiltonian exactly within a finite size cubic supercell 
with periodic boundary conditions. In the supercell, Mn 
and As anti-sites are randomly distributed within the 
lattice and Mn spins are described by the classical 5/2 


FIG. 1. (Color online) Wave function probability distribution 
of a single hole with a single Mn impurity in our supercell 
(black cubic outline). The lattice structure shown depicts the 
actual lattice constant with respect to the size of the supercell. 
The light color sphere in the center of the cloud is the Mn 
impurity. 


III. RESULTS AND DISCUSSIONS 

A. Low doping limit: 
bound holes &; pair-bonding states 

The advantage of our method is the ability to cap¬ 
ture the low Mn doping as well as high Mn doping limit. 
In this section, we first study the low limit where there 
is only a single Mn impurity in our simulation super¬ 
cell (0.02% Mn). Figure [l] clearly shows a hydrogenic- 
like bound state with a Bohr radius of ^ 4 aj c , where 
a ic = 0.565 nm is the lattice constant of GaAs. The 
outer box shows our simulation supercell, and the lattice 
structure inside follows the actual lattice constant with 
respect to the box. The position of the Mn, shown as 
the light sphere in the middle of the distribution, is ran¬ 
domly chosen. The binding energy of this bound state is 
^ 40 meV, which is about 1/3 of the experimental value 
for isolated Mn. The results approach the experimental 
value with an increase in size of our simulation cell and 
the number of k points for the envelope function expan¬ 
sion. 

If two bound states are placed close to one another, 
one expects that bonding and anti-bonding states will be 
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FIG. 2. (Color online) Wave function probability distribution for quasiparticle state with two Mn impurities sitting in the first 
to fourth nearest neighbor with respect to each other. The outer boxes are not the exact size of our simulation super cells. 
Light color spheres are the Mn and the small inner cubics show the actual size of unit cell of GaAs. (a) First nearest neighbor 
Mn with binding energy 276 meV, (b) Second nearest neighbor Mn with binding energy 232 meV, (c) Third nearest neighbor 
Mn with binding energy 222 meV, (d) Fourth nearest neighbor Mn with binding energy 199 meV. 


formed due to the interaction between the states. De¬ 
pending on the strength of interaction due to this over¬ 
lap, the energy difference can be very large. In our case, 
the bonding state will move deeper in the gap while the 
anti-bonding state will move to the valance band region. 
Figure [2] shows our simulation results for the bonding 
state of two Mn impurities. Arranged from first to fourth 
nearest neighbor positions, the binding energies are 276 
meV, 232 meV, 222 meV, and 199 meV respectively. To 
provide a better visualization, the images show a zoomed- 
in view of the probability distribution function, i.e., the 
outer box outline is not the actual size of our simulation; 
however, the inner cubic outline is the actual size of the 
lattice structure with respect to the probability distribu¬ 
tion. The highly anisotropic structure of the probability 
distribution is consistent with STM experiments.^ 3 As 
expected, the bonding states have much larger binding 
energies than the isolated bound states. 


B. High doping regime 

1. DoS and LDoS 

The impurity band model assumes that a detached 
band still persists in the high doping limit. In this sec¬ 
tion, we evaluate this high doping regime by studying 
the density of states (DoS) and local density of states 
(LDoS) in order to confirm whether or not the detached 
impurity band persists. As we know from the previous 
discussion, two nearby bound states will form bonding 
and anti-bonding states. In the high doping limit, there 
are many impurity pairs which give rise to a great amount 
of bonding states in the gap region. On the other hand, 
the anti-bonding states will contribute a large number of 
states to the valance band region. 

The total DoS can capture the distribution of these 
states. Figure |3] shows the DoS plots for 5% Mn doping, 


without As anti-site compensation. This figure shows 
the DoS with and without exchange interactions, i.e. 
ferromagnetic and non-ferromagnetic phases, as well as 
the DoS obtained from virtual crystal approximation 
(VCA). 26 VC A assumes the wave vectors k remain good 
quantum numbers with disorder treated as an energy 
spectrum broadening. As shown in the DoS plot, the 
VCA can only change the structure of the valance band 
by slightly shifting the energy levels into the gap, whereas 
our more realistic treatment of disorder tends to pull 
states near the top of valance band deeper into the gap 
region. The kinetic-exchange interaction further splits 
energy levels and redistributes the states in the gap. As 
expected, there is an extensive weight of the DoS in the 
mid-gap region which arises from the deep bonding states 
created by neighboring impurities. However, these DoS 
curves show only a disordered valance band as opposed 
to an impurity band detached from the top of the valence 
band. 

To further study the spatial inhomogeneity of the dis¬ 
tribution of states in the gap region, we calculate the 
LDoS which can be qualitatively compared with the STM 
experiments. Figure [4] shows our simulation results along 
a line in our simulation cell. The gray vertical plane is the 
Fermi energy. The peaks formed by bonding states in the 
gap region are consistent with the STM measumentsj^ 
which indicates our model captures the disorder effects 
correctly. We emphasize here that in our approximation, 
which ignores the Fock-exchange contribution, we cannot 
capture the expected reduction of the DoS at the Fermi 
energy in disorder systems. 


2. Localization properties of the states in the high doping 
regime 

From our DoS study, an impurity band detached from 
the valence band does not persist for high Mn doping. To 
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FIG. 3. (Color online) Density of States (DoS) with and with¬ 
out exchange interactions for 5% Mn doping without As com¬ 
pensations. We shift all the curves to make the Fermi ener¬ 
gies to be zero as the vertical line shows. The red curve is 
the DoS obtained from VC A with an energy level broaden¬ 
ing of 60 meV. The inset shows one typical Mn distribution 
in our simulation cell. The purple line in the cell shows the 
coordinate of the LDoS calculation (shown in Fig. 4). 



FIG. 4. (Color online) Local Density of States (LDoS) for 
5% Mn density without As compensations. This plot shows 
different positions along the x-axis of our simulation supercell 
which is shown in Fig. [ 3 ] The vertical plane shows the posi¬ 
tion of the Fermi energy. The results agree qualitatively with 
the experiments in Ref. IT2 , 

further confirm this observation, we study the localiza¬ 
tion character of the states throughout the spectrum. If 
an impurity band were to exist in the high doping limit, 
we would expect localized states in both tails of the im¬ 
purity band. As we have shown, bound states in the low 
Mn doping limit (< 0.1%) are well localized indicating 
the existence of an impurity band for this limit. The 
probability distribution of these wave functions are con¬ 
centrated in a small region with a size of several GaAs 
lattice constants. When increasing Mn doping, some of 
the hole states which are well isolated from others will 
stay localized, but for some others their wave functions 
begin to overlap with each other leading to delocaliza¬ 
tion. To have a quantitative description of localized and 


delocalized states, we use participation ratios (PRs) to 
characterize the extent of localization of states and com¬ 
pare with the probability distributions directly to estab¬ 
lish the criteria for localized states. Then, we group these 
states as localized and delocalized based on this criteria. 
The results show that our model is able to capture the 
localization physics with its more realistic treatment of 
disorder. PR is defined as the inverse of the integral over 
the si mulat ion volume of the wave function to the fourth 

power 


PR vf v d^mf)r (2) 

where V is the volume of the simulation cell. T(r) is the 
normalized wave function. A simple example gives the 
meaning of the PR. For delocalized states, |\h(r) | 2 ^ l/W, 
so PR ~ 1, but for localized states, the wave function 
probability is concentrated in an area much smaller than 
V. The PR has the order of 1/F, which is much smaller 
than 1 if V is large enough. Another way of distinguishing 
localized and delocalized states is using size scaling of the 
PRs. With this method, the PR is almost constant for 
delocalized states, but scales as 1/V for localized states. 
The disadvantage of using size scaling of PR is the com¬ 
putation time. In this paper, similarly to other studies 
in disorder systems, we use the magnitude of PR and re¬ 
late it to the probability distribution for hole state wave 
functions to distinguish localized and delocalized states. 

Figure [5] shows the participation ratios as a function of 
energy and the typical wave function probability distri¬ 
butions with different participation ratios. For our choice 
of V, states with PR <0.1 are localized states with their 
wave function concentrated in a small area. For PR >0.1 
but close to 0.1, these states are a transition from local¬ 
ized to delocalized states which define the mobility edge. 
Within this definition there is of course no sharp sep¬ 
aration of localized and delocalized states. This would 
require the use of the finite-size scaling method. How¬ 
ever, for the purposes of our study it suffices to treat 
states with PR< 0.1 as localized and above it as delocal¬ 
ized. With this criteria, the red region in Fig. [6] shows 
the area of the localized states in the DoS plot, mean¬ 
ing that states deep in the gap are well localized which 
is consistent with our previous discussion. These states 
are bonding states formed by nearby Mn impurity pairs 
with strong interactions. However, if more Mn are doped 
into this system, localized holes have more places to hop 
around becoming delocalized due to the increasing num¬ 
ber of nearby Mn. There will be less localized states for 
higher Mn doping samples. Figure [7] shows the simulated 
results for different Mn doping without compensation ob¬ 
tained by averaging over four disorder realizations. It is 
clear that the number of localized states decreases with 
increasing Mn doping. In our results only states deep in 
the gap region are localized far away from the Fermi en¬ 
ergy, which does not support the impurity band picture. 
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FIG. 5. (Color online) (a) Wave function probability distribution with PR = 0.015, (b) Wave function probability distribu¬ 
tion with PR = 0.1, (c) Wave function probability distribution with PR = 0.2, (green spheres denotes Mn impurities), (d) 
Participation ratios. 
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FIG. 6. (Color online) Density of States (DoS) for 5% Mn 
without As compensations. Red area shows the localized 
states 


C. Ac-conductivity 

In this part, we will study the ac-conductivity to ex¬ 
plain the spectroscopic experimental observation which 
was used to support the impurity band model. The 
method used here is the standard linear response Kubo 
formalism. Early experiments of infrared ac-conductivity 
measurements on Gai_ x Mn x As shows a broad peak near 
200 meV for metallic samples. A red shift was ob¬ 
served when increasing hole concentration, which is the 
key evidence of impurity band model. 7 Later on, similar 
experiments on optimized growth and post-growth an¬ 
nealing samples show a blue shift when increasing hole 
concentration.^ By post-growth annealing procedures, 
the density of compensating defects and other uninten¬ 
tional impurities can be greatly reduced. To compare 
with these experiments, we will study uncompensated 
samples first. Figure [8] shows the ac-conductivity for Mn 
dopings ranging from 1% to 5% without compensation. 
For each doping rate, we do see a broad peak in the low 
energy region ^100 meV. When increasing Mn doping, 
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FIG. 7. Average number of localized states in our simulation 
cell. Data obtained by averaging over 4 disorder realizations. 

there is a red shift when Mn is less than 3% and no shift 
for higher dopings. We expect that compensation plays 
an important role. Aside from the peak ~100 meV, there 
is one other peak ^800 meV for lower Mn doping sam¬ 
ples, which arises from an unphysical finite size effect and 
band cut-off for the purposes of the calculation. To study 
the effects of compensation, we fix the Mn doping to be 
5% and consider As anti-sites as a source of hole compen¬ 
sation. Figure [9] shows our simulated results for different 
hole densities. Compared to the uncompensated samples, 
the peaks are now shifted from ~100 meV to ^200 meV. 
There is still no prevailing red shift with increasing hole 
density. 



Energy [eV] 

FIG. 8. (Color online)Real part of a xx for various Mn density 
without As compensations. Data obtained by averaging over 
4 disorder realizations. 

The simulated results show our model is able to cap¬ 
ture the mid-infrared peak in optical absorption mea¬ 
surements by considering the transitions from states be¬ 
low the Fermi energy to states above the Fermi energy 
and deep in the gap region. The Fermi energy resides 
in the delocalized disordered GaAs valance band region. 
This demonstrates that the assumption of the detached 


IB model is not valid from its starting point in such 
metallic samples. As shown in the ac-conductivity simu¬ 
lation, the peaks for compensated samples move towards 
high energy compared to uncompensated samples. Hence 
one can get any pattern of the peak shifts if less and 
more compensated (annealed and unannealed) samples 
are mixed. Experimentally, the actual density of com¬ 
pensation defects is hard to determine. Aside from As 
anti-sites, there are other types of charge and moment 
compensations. So, based on our simulations, the peak 
experimentally observed shifts of the ^200 meV cannot 
be used to make any direct conclusions regarding the na¬ 
ture of states at the Fermi energy. 
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FIG. 9. (Color online) Real part of a xx for 5% Mn doping 
and different hole density due to As compensation. Data is 
obtained by averaging over 4 disorder realizations. 


IV. CONCLUSIONS 

In this paper, we provide a numerical study of the 
effect of disorder and localization on the states near 
the top of the valence band and in the gap region of 
Mn-doped GaAs. Through exact diagonalization of our 
model Hamiltonian, the low and high doping limits are 
captured. The local spectroscopic properties are stud¬ 
ied and compared with STM experiments by simulating 
the DoS and LDoS. Furthermore, via visualization of the 
probability distribution of the quasiparticles and the con¬ 
nection to their respective participation ratios, we sepa¬ 
rate the states into localized and delocalized. It is shown 
for metallic samples that the states deep in the energy 
gap are well localized and states near the Fermi energy 
are delocalized. The optical conductivity with different 
Mn dopings and As compensations are then simulated. 
Our results show that the peak at ^200 meV is due to the 
transition of states below Fermi energy to states in the 
gap region. There is no direct relation between the opti¬ 
cal absorption measurements and the nature of states at 
the Fermi energy. Our overall results are consistent with 
the disordered valence band model and disagree with the 
assumption of the IB models of a detached impurity band 
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from the valence band. 
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